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§1. SOLITONS AND THE SCHRODINGER OPERATOR 

The modern theory of exactly solvable one-dimensional and two-dimensional 
Schrodinger operators 

L=-d 2 x +u{x) (n = l), 
L=(B + B(x,y))(d + A{x,y))+2V(x,y), d = d x - id y , z = x + iy (n = 2) 

cannot be separated from the theory of spatial one- and two-dimensional integrable 
nonlinear systems of soliton theory. This is especially so in the rapidly decaying 
and periodic case (when the coefficients of the operators are rapidly decaying and 
periodic functions of the space variables) . 

I. One-dimensional case (n = 1). In the one-dimensional case, it follows from 
the papers [GGKM, L], where the famous KdV equation 

u t = 6uu x — u xxx 

was solved in the rapidly decaying case that it can be viewed as a kind of "sym- 
metry" of the spectral theory of the class of operators L = —d x + u(x) acting in 
the Hilbert space £2 (R) • This means that the KdV equation can be written in the 
form 

L t = AL — LA = [A, L], 

where A is the linear differential operator A = —Ad^. + 3(ud x + d x u) . Therefore, for 
example, the eigenvalues Lip = Xiip of the operator L on any translation-invariant 
class of functions of x on the line R turn out to be integrals of the KdV equa- 
tion: dXi /dt — . Using the previously solved inverse scattering problem, KdV was 
solved in [GGKM] for rapidly decaying functions. Beginning with S. P. Novikov's 
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paper [N74], in the articles [DN74, D75, IM, L75, MV] the solution of the peri- 
odic problems for nonlinear KdV from soliton theory was obtained simultaneously 
with that of the inverse problems for the spectral theory of the linear operator 
L , which had not been previously solved in the periodic case (see the surveys 
[DMN, Kr77, DKN2] and the book [Sol]). The higher analogs of KdV yield an 
infinite-dimensional commutative group of such "spectral symmetries" for the op- 
erator L . The finite-dimensional orbits of this group generate the famous rapidly 
decaying and quasiperiodic ( "multi-soliton and finite-zone") solutions of KdV as 
well as the remarkable exactly solvable Schrodinger operators on the line in £ 2 (R) ■ 
Here the hyperelliptic -functions of Riemann surfaces make their appearance, to- 
gether with a series of other now widely known formulas and results. Such are the 
consequences of continuous spectral symmetry. 

A few years ago some interesting advances were made [W, SY, SI] in the applica- 
tion of discrete symmetries, such as sequences of the well-known Backlund-Darboux 
transformations, to the spectral theory of the operator L = —d% + u(x) . Probably 
the most significant were obtained by A. B. Shabat and A. P. Veselov [VS1]. 

II. Two-dimensional topologically trivial case (n = 2). Beginning with the 
papers [M, DKN1] in 1976, as the basis of the two-dimensional analog of the theory 
of integrable systems associated with the two-dimensional Schrodinger operator, it 
is customary to take a representation of the form 

L t = LA + BL, 

where L is the two-dimensional Schrodinger operator (see above), A and B are 
differential operators on the plane (x, y) . Here the data of the inverse problem for 
the spectral theory of the linear operator L are taken from one level only, namely 
from the set of solutions of Lip = under some "boundary" condition or other 
on ip . For the operator L with doubly-periodic coefficients one takes the Bloch 
waves 

T a ip = exp{ip a T a }ip 7 Lip = 0, a = 1, 2. 

Here T a are the periods, a = 1,2, while T a are the translation operators by the 
period T a , and p a is the quasi-momentum. The set of all complex Bloch waves of 
zero level Lip = constitute a one-dimensional complex manifold T: ip(P,z,z), 
P G r (the Fermi curve). Without going into the details of the analytic properties of 
the function ip , first clarified in [DKN1], let us note the following: exactly solvable 
topologically trivial two-dimensional Schrodinger operators correspond to the case 
in which one Fermi curve T is of finite genus. The selection of the inverse problem 
data corresponding to the purely potential operators 

L = dd + V(x,y) 

was carried out in [NV1] (see the survey [NV2]). The simplest two-dimensional 
system of soliton theory, generalizing KdV, is of the form 

^ = (d 3 + d 3 )V + d(uV) +d(uV), du = 3dV 
and can be represented as 

BT ^ 
Tt =%B] +f L, 

B = (d 3 + d 3 ) + (ud + ud), f^du + du. 
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There is an entire infinite-dimensional commutative group of systems of this form 
(the higher analogs of two-dimensional KdV or Novikov-Veselov hierarchy). 

Such is the two-dimensional analog of the continuous spectral symmetry of the 
class of Schrodinger operators. In contrast to the one-dimensional case, it acts only 
on the eigenfunctions of one level Lip = . 

Already at the end of 18th century, Laplace (see [L]) proposed substitutions 
acting on the solutions of the equation Lip — and transforming them into solutions 
of a new Schrodinger equation Lip — . We call them "Laplace transformations" . 
In the present paper we use the Laplace transformations as a kind of "discrete 
symmetry" for the spectral theory of the Schrodinger operator L corresponding to 
one level Lip = , i.e., being in good agreement with the ideology of two-dimensional 
soliton theory described above. 

More precisely, we shall construct new classes of two-dimensional Schrodinger 
operators that are exactly solvable at one or several energy levels also in the topo- 
logically nontrivial case, a case about which little is known (practically nothing, 
except the papers [AC, DN80] and the survey of one of the authors [N83]). It is 
precisely this last part of the present paper that we regard as its main result. Part 
of the results of the present paper was already announced in [NV3] . 

An alternative approach, based not only on Laplace trancformations, but also on 
additional symmetries, is being developed by A.B.Shabat and A.P.Veselov in [S2] , 
[SV2] (in preparation). Its relation with spectral theory of topologically nontrivial 
operators (see below) is totally unclear. 

§2. Two-dimensional Schrodinger operator in a magnetic field, 
topologically nontrivial operators. the pauli operator 

The general two-dimensional Schrodinger operator on the Euclidean plane R 2 
has the form 

L = (8 + B)(d + A) + 2V = (d + A)(B + B) + 2U, (1) 

where A, B , V , and U are functions of (x, y) , z = x + iy . It is characterized by 
the electric potential V(x, y) and the magnetic field (directed along the third axis): 

2H(x, y)=B z -A Y , U = V + H (2) 

up to gauge transformations 

L — > e*Le~, ip — > e' ip, 

where f(x,y) is a smooth function and ip G /^(R 2 ) • 

We shall consider the periodic case, when the group of periods is given on the 
plane, this group being generated by the basis vectors T i and T 2 , so that the 
magnetic field and the electric potential are doubly periodic: 

H(x+T 1 ) = H(x) = H(x+T 2 ), x = (x,y), 
V(S+%) = V(x) = 0V(x+T 2 ). 

Denote by K the elementary cell in the (x, y) -plane, i.e., the interior of the 
parallelogram with vertices 

K = {0,0+T 1 ,0+T 1 +T 2 }. 
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By the magnetic flux we mean the expression: 

[H] = JJ^H(x,y)dxdy. 

Definitions, a) The case in which [H] = is called topologically trivial. Here all 
the coefficients of the operator L are doubly-periodic. 

b) The case in which [H] =/= is called topologically nontrivial. The coefficients 
of the operator L are not doubly-periodic functions on the (a;, y) -plane. Special 
attention will be paid to the case in which [H] = 2im , where n is an integer. Wc 
shall say that this is the integer case. 

Suppose, for the sake of simplicity, that the lattice is rectangular, 

I*i = (71,0), T 2 = (0,T 2 ). 

Definition. The magnetic translations are the following transformations: 

Tiip(x,y) = exp{if 1 }ip(x + Ti,y), 

T 2 tp(x,y) = exp{if 2 }tp(x + T 2 ,y), 

where 

A(x + T 1 ,y)-A(x,y)=idf 1 , B(x + T 1 ,y) - B(x,y) 
A(x, y + T 2 )- A(x, y) = idfc, B{x, y + T 2 )- B(x, y) 

It is known that magnetic translations commute with the operators L, i.e., 
LTi — TiL. They also satisfy the following commutation relation: 

f 2 fi = e^^fifa, (4) 

where [H] is the magnetic flux through the elementary cell K . These properties 
can be verified by an elementary substitution. 

Thus relation (4) implies that only in the integer case [H] — 2nn can we define 
the magnetic Block functions: 

T k ip(x,p) = exp(ip k T k )ip(x,p), k = 1,2. (5) 

Here the magnetic quasi- momenta (pi,p 2 ) range over the points of the two- 
dimensional torus p e T"l corresponding to the reciprocal lattice . For fixed 
quasi- momentum (pi,P2), the operator L acts as an elliptic operator on the sec- 
tions of the complex linear bundle over the torus with the connection whose 
curvature H{x,y) is determined by formula (5). If the magnetic field, the electric 
potential, and the quasi- momenta (pi,p 2 ) are real, then this operator is semi- 
bounded, selfadjoint, and has a discrete spectrum of finite multiplicity 

Lipj = -Sjipj, Sj = e j (p 1 ,p 2 ), j = 0, 1,2, . . . , s < £\ ^ s 2 < . . . , 

for given (pi,p 2 ) . If the energy is large (e — > oo), then the sizes of the magnetic 
zones 

A,- = max £.,■(») — min £,(») 

tend to zero: Aj — > as j — > oo . It is assumed that the magnetic flux is nonzero: 
[H]^Q. 

It is of interest to study the topology of generic operators, initiated in [N81] (also 
see the survey [N83]). 



(3) 



idfi, 
idf 2 . 
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Definition. A selfadjoint operator L with smooth doubly-periodic (V, H) , [H] = 
2im is called generic if for all j and all (pi,P2) € T 2 we have 

£j{P\,P2) + £k(pi,P2), j^k. 

In this case the eigenfunctions Lipj = £jipj in the Hilbcrt space Ti of sections 
of the bundle over the (a;, y) -torus T 2 depend on the quasi- moment a (pi,P2) € 
as on parameters. This generates a linear complex bundle 77 j over the torus 
T 2 , depending on the number j . The corresponding Chern classes c\ (rjj ) , whose 
study was initiated in [N81, N83], as was later established by physicists [T], play 
a fundamental role in the theory of the integer quantum Hall effect, discovered in 
1981. These classes may assume different random series of values. 

Problem. Study the asymptotic properties of the integers C\{r)j) as j — > oo for 
generic operators L when [H] = 2nn . 

The eigenfunction ipj(x,y,pi,p2) , where Lipj = £j(pi,P2)ipj , in the generic case 
possesses a manifold of zeros 

ipj(x,V,Pi,P2) = 0. 

This equation expresses two real conditions. Hence, generally speaking, the man- 
ifold of zeros Mj is a two-dimensional real submanifold in the Cartesian product 
of tori 

Mj c T 2 x T 2 , Mj = Mj U Mj, 

where Mj is projected into a finite number of points of the torus T 2 , while Mj at 
a generic point is given locally in the form of a graph 

p 1 =F jl (x,y), p2 = F j2 (x,y). 

The problem of studying the analytic properties of the collection of vector func- 
tions p = Fj(x,y) was posed in [N83]. It seems probable that the problem of the 
spectral theory of the operator L can be effectively solved on the basis of these 
notions. 

Of special interest is the nonrelativistic Pauli operator for spin 1 /2 on the plane 
with an electric field parallel to the plane and a magnetic field perpendicular to it. 
This is a vector operator. It has the form (the charge, mass, and Plank constant 
are assumed equal to 1 ) 

2 

P = Y. {d *-iA a ) 2 + (W + a 3 H), (6) 

where x\ = x, x 2 — y, d a — d/d a , while (Ai,A 2 ) is the vector potential, W is 
the electric potential, H = d\A 2 — d 2 A\ is the magnetic field, and 
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The operator P has the form of a direct sum of two scalar operators 

P = L+ L-, 

2 2 

L + = Y J {d a -iA a ) 2 + {W + H) 1 L_ = Y,(d*-iA a ) 2 + (W-H). 

a—l a—1 

For a zero electric field W = , the Pauli operator may be factorized as 

L+ = (d + A)(B + B) = (d + B)(d + A) + 2H, 

L_ = (d + B){8 + A) = (d + A)(d + B) - 2H. 

It is convenient to choose a "real" gauge B = —A, A = A 2 + iA x . We shall also 
impose the usual Lorentz condition d\Ai + diAi = or 

Im(dA) = Q. (9) 

According to the results of [AC, DN80], in the absence of any electric field, the 
Pauli operator (-P), which is always nonnegative under these conditions, has the 
point e = as its ground state for rapidly decaying fields H (see [AC] ) and doubly- 
periodic H when [H] ^ (see [DN80]). Here if [H] > 0, then the ground state 
corresponds to the + sector and satisfies the first order equation 

{d + B)ip = 0. 

The equation L-ip = in this case has no solutions with the necessary boundary 
conditions. If [H] = 2im, n > 0, then the basis of the magnetic Bloch states may 
be written in the form 

n 

*P = e' P Y[o-(z-a j )e az , aj eK, 

where 

A(p = —H, ip = — J I In \a(z — z')\H(z') dxdy, z = x + iy; 
kJJk 

here a\, . . . ,a n are arbitrary constants, while a is determined by them from the 
condition requiring the quasi- momenta to be real (see [DN80]). 

The degree of degeneracy for the level e = is the same as that for the Landau 
level for H = const . The magnetic field enters into the formula only via <p and a . 

§3. Laplace transformations. Main results 
(the topologically trivial case) 

For the Schrodinger operator L = (d+B)(d+A) + 2V , we start with an arbitrary 
solution of Lip = . Then the function = (d + A) tp satisfies the equation 

[V(d + A)V- 1 (B + B) + 2V] $ = 0. 

Thus we see that if Lip = , then Lip = , where 

L = V(d + A)V- 1 (d + B) + 2V, ^={d + A)ip. 

Since Vid + A^V- 1 = d + A, where A = A-(lnV) z , we have L = (d + B){d + A) + 
2V, where B = B, A = A — (In V) z , 

[H = H+\dd{\nV), 
[V = V + H. 

It is essential that the final formulas (10) provide the transformations only in 
terms of physical quantities. 
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Lemma 1. Suppose that the potential V(x,y) and the magnetic field H(x,y) are 
real and doubly-periodic, and 7^0. In this case the magnetic fluxes [H] and [H] 
coincide, while for the potentials we have the relation 

[V] = [V\ + [H\. 

The proof follows from formula (10) for the Laplace transformations because 
[dd(lnV)} = 0. 

Lemma 2. Let the potential V and the magnetic field H be real, while the operator 
L is given in the real Lorentz gauge 

L = (d + B)(d + A) + 2V, B = -A, Im^- = 0. 

Then for any solution of Lip — the function ip — e~®(d + A)ip satisfies the 

equation Lip — 0, where L — (d + B)(d + A) + 2V , and the operator L is gauge 
equivalent to the image of the Laplace transformation 

L = e- Q Le +Q , V = V, B = B + Q-, A = A + Q Z . 

Now if Q = In \fV , then the new operator L can be expressed in the real Lorentz 
gauge as follows 

B = B + \{\nV)z, A = A-\{\nV) z , —B = A, ImA* = 0. (11) 

The proof of this lemma also follows from the formulas specifying the Laplace 
transformations. 

Speaking of the Laplace transformation for magnetic fields and potentials, we 
shall not distinguish, as a rule, operators from the same gauge equivalence class. 
But if we are concerned specifically with real selfadjoint operators, then by Laplace 

transformations we mean the transformations L — > L effected in accordance to 
formulas (11). 

Now let us consider infinite sequences of Laplace transformations: suppose that 
Vj = cxp fj , then 

e /3+1 - e f i + H j+1 , H j+1 = H, + \ddf, , 

whence 

\ Afj = e f i +1 - 2e^ + e 1 ^ 1 . (12) 

Lemma 3. After the substitution fj = ipj — tpj-i , the infinite sequence of Laplace 
transformations reduces to the well-known two-dimensional Toda lattice 

|A^- = e^+ 1- ^ - eww- 1 . (13) 

Indeed, the substitution fj = ipj — ifj-i in (12) yields 

iA<^ = e w+i-w - e <Pi-<ei-->- + h{z,z), 
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where h(z,z) is a function that does not depend on the number j . It is easy to 
get rid of it by using the change ipj — ► ipj + a(z,z) , where (1/2) Aa = h. 

The two-dimensional Toda lattice together with the representation in the form of 
an (L-A)-pair was found by A. Mikhailov in [Mik] and, independently, in [LS, BUL] 
for finite chains with free ends, as a generalization of the Liouville equation related 
to Lie algebras. It was known to classical geometers (see [Dar, Tzi]) as a sequence 
of Laplace transforms. The connection of classical geometric studies with systems 
from soliton theory has been observed by A. M. Vassiliev (unpublished), but no 
attempts to use this observation were made. A consequence of this comparison is 
the following 

Proposition 1. The algebro- geometric solutions (expressed in terms of 6 -functions)^ 
found in [Kr81] for the two-dimensional Toda lattice generate for any j a two- 
dimensional Schrodinger operator Lj with quasiperiodic coefficients, which pos- 
sesses a Bloch solution Ljtp = with appropriate analytic properties on the same 
Riemann surface Y of finite genus. Conversely, suppose that we are given a non- 
singular algebraic curve Y on which two distinguished points P + , P_ with local 
parameters w + , W- are fixed, together with the divisor D = Pi + ■ ■ ■ + P g of 
degree g . If the operator L = Lq corresponds to the data of the inverse problem 
(r, P±, w±, Pi, . . . , P g ) , then the operators Lj from the chain of Laplace trans- 
formations are obtained from the data (T, P±, w±, P{, . . . , Pj>) , and we have the 
following linear equivalence of divisors 

Pi + ■ ■ ■ + PI ~ j(P+ - P-) + Pl + ■ ■ ■ + Pg- (14) 

Recall (see [DKN] ) that in the generic case P\ , . . . , P g are poles of order one of 
the Bloch solution of Lip = that do not depend on x, y. Near the points P± 
the function ip has the asymptotics 

V>~ Cl (x,y) e z ^+(l +a(w+)), P-P+, 

z = x + iy. (15) 

4>~c 2 (x,y)e z / w -(l +*(w-)), P~P_, 

The function is meromorphic on Y \ (P + U P_) and depends on x, y as on 
parameters. Such a function always satisfies an equation of the form Lip = 0, 
where 

L = (d + B)(d + A) + 2V, A=-d\nc 2 , B = -dlnci. 

Indeed, since (d + A) c 2 = (d + B) c\ = , the function (d + B)(d + A)tp possesses 
similar analytic properties and, therefore, differs from ip only by a constant factor. 
It is easy to see that the function ip = (d + A) ip has the same poles Pi , . . . , P g and 
asymptotics in P± of the form 

ip^c 1 (x,y)w~ 1 e z /^(l + a( W+ )), P^P+, 
c 2 (x 1 y)w-e z/w -(l +a(w-)), P~P_. 



More generally, we have the following 
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Lemma 4. The function i})^ > being the result of j applications of the Laplace 
transformation to the operator L , besides its poles P\ , . . . , P g , possesses singulari- 
ties at the points P± of the form 

^ =c[ l) (x 1 y)w^e z ^(l + a(w + )), P ~ P +1 
^0) = c ^(x,y)w j e z / w -(l + o-(w-)), P~P_. 

The proof of Proposition 1 now follows if we compare the result of Lemma 4 
with the construction of the algebraic solutions of the two-dimensional Toda lattice 
proposed by I. M. Krichever [Kr81]. 

Explicit formulas involving O -functions for the coefficients of the corresponding 
operators Lj can be found in [D81, p. 50] and [Kr81, p. 77]. They imply, in 
particular, that the coefficients are quasi-periodic functions of x and y , so that the 
operators Lj are topologically trivial. 

The solutions of a periodic chain of period N , <p n +N = , correspond to curves 
with points P + , P_ such that NP + — NP- ~ 0, i.e., curves on which there exists 
a meromorphic function with a unique pole of iVth order at P+ and a zero of order 
N at P_ . To such solutions correspond cyclic chains of Laplace transformations. 
It turns out that under certain analytical assumptions all such chains are algebro- 
geometric, namely, we have the following result. 

Proposition 2. Suppose that a chain of Laplace transformations is cyclic, V^+jv = 
Vi , and all the potentials Vi — exp /j are doubly periodic, smooth, real, and positive 
functions. Then for all n the Schrodinger operators L n are algebro- geometric with 
respect to the zero energy level L n ^ = . 

For N = 1 from the relation (12) we get A/i = , which in the doubly periodic 
case implies f\ — const . The case N = 2 reduces to a similar question for the 
sinh -Gordon equation 

Aip + sinhtp = 0. 

The corresponding result in this case was first obtained by Pinkall and Sterling 
[PS]. A. Bobenko considerably simplified the proof [B]. The proof of Proposition 1 
presented below is a direct generalization of Bobenko's arguments. 

It follows from Lemma 3 that cyclic chains of Laplace transformations correspond 
to doubly-periodic solutions of the periodic Toda lattices (13). It is known (see 
[Mik]) that such a chain may be represented as an (L- A) -pair of the form 

[d + P,d + Q] = 0, (18) 

where P and Q are the following matrices depending on the "spectral" parame- 
ter A: 

A-^-^X 


e ¥>n-l-¥>n / 

This allows to determine, in the standard way, an infinite series of commuting 
flows 

[d + P,d i + P i ] = 0=\B + Q,d i + P i \, (19) 



P = 



(dip! 1 


V A 



\ 

1 



Q 



( 



V o 
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where di = d/dU and Pi is a matrix polynomially depending on A (see, for 
example, [DS]). The derivatives £W = (d/dti)(ip) satisfy for all i the linearized 
system 

lA^ l) = £, { h i l 1 e ipk + 1 - Vk -^(e^+i-Vfe + e ^-^-i) + ^W^^-Vfe-i. ( 20 ) 

This system is an elliptic linear system on a compact manifold (the torus T 2 ) , 
hence there exists only a finite number of linearly independent £W . This in its 
turn implies that any doubly-periodic solutions of the Toda lattice is stationary 
for an appropriate higher flow (19) and by a standard argument (see [DKN2]) is 
algebro-geometric. This concludes the proof of Proposition 2. 

§4. Laplace transformations and 
topologically nontrivial periodic operators 

As indicated above, in the nonsingular periodic case, the Laplace transformation 

H k+1 =H k + iA/ fe , e'* = Vfc, e^ 1 = eJ« + H k+1 

preserves the magnetic flow [Hk+i] — [H k ] and changes the mean value of the 
potential 

[ e /*+i] = [ e A] + [H k+1 ]. 
Thus (nonsingular) cyclic chains with nonzero magnetic flux do not exist. 
Definition. A semi-cyclic chain of length n is a chain such that 

H n = H 0l V n = V + C n . (21) 

In the nonsingular case we obviously have 

C n =n[H }. (22) 

Here we assume that all the operators in the chain L , L\, . . . , L n have nonsingular 
coefficients and the invariant expressions are doubly-periodic 

H n {x + Ti) = H n (x + T 2 ) = H n (x), 
V n (x + T 1 ) = V n (x + T 2 ) = V n (x). 

Here T\ , T 2 are the basis vectors of the lattice of periods in the plane, and x — 
(x 1 ,x 2 ). 

Example 1. Let n = 1. From the formulas for the Laplace transformations, we 
immediately obtain 

H ! = #o, e h = e fo + ffi , A/ = 0. (23) 
In the nonsingular doubly-periodic case, we have 

/i = const, /o = const, H = const . 
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Thus the operator L\ is of the form 

Li=L + Ci. (24) 

Here L\ is the Landau operator in a homogeneous magnetic field and trivial (con- 
stant) potential. If we know solutions of ioV'o = in the Hilbert space £ 2 (R 2 ) , 
then the function ipi = Qoipa = e ta / 2 {d + A a )ip gives the solution ipi £ £ 2 (R 2 ) • 
Iterating this procedure, we get 

L n = L + nC u = Qn-i ■ ■ ■ Qo{ifo) e C 2 (M. 2 ). (25) 

In this case we obtain the eigenfunctions of the operator L with eigenvalues 
(-nCi) in the space £ 2 (M 2 ) 

£o^n = (-nCi)^ n , n = 0,1,2,.... (26) 

It is easy to see that 

d = H = [H ]/K = const, (27) 

where K is the volume of the elementary cell. We know that the Schrodingcr 
operator (— Lq) is semi-bounded. There exists a constant E such that 

{-L^^^-Eo^,^), ^e£ 2 (M 2 ). 

Hence the solution LoV'o = belonging to £ 2 (IR 2 ) is possible only in the case 
[H ] > 0. If W = (or V = H ), i.e., if 

L = {d + A)(B + B) = (d x + Mi) 2 + (d y + iA 2 f, 

then the solutions serving £ 2 (R 2 ) can be found in the form (d + B)tpo — for the 
real Lorentz gauge B = — A, A\ x + A 2y = 0. (This will be established below in 
connection with the theory of quasi-cyclic chains.) In this case the spectrum of the 
operator (— Lq) is discrete and has the levels 

A„ = nCi = nH , (-L )ip n = \ n ip n . 

This is precisely the spectrum of the Landau operator (up to a translation by the 
constant H /2). However, if [H ] < 0, then we must begin with the operator 

L' = (d + B)(d + A), V = 

and search for the solution in the form (d + A) tf>' = . The spectrum is obtained 
by applying the inverse Laplace transformation (and its iterations): 

Vi=QoV>o = e W2 (9 + A)^, go^lnVKo, 
L\ =L' + C[, C[ = -H X n = -nH , {-L' )^ n = A„Vn- 

This concludes Example 1. 

Our main results are related to the following class of operators. 
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Definition. A quasi-cyclic chain of Laplace transformations of length n is a chain 
Lo, Li, . . . , L n such that 

C n + H n = V n = e f " , H = V = e fo . (28) 



Suppose further that the magnetic flux is positive and integer- valued, i.e., [H ] = 
2irm > . 

Example 2. Let n = 2. Let us write out the equations that follow from the 
conditions of semi-cyclicity and quasi-cyclicity: 

^Ag = — C 2 — 4e a / 2 sinhg, g = fo — a/2 (semi-cyclic case), (29) 
|A/ = C 2 - 2e /o , e h =2C 2 -e fo (quasi-cyclic case). (30) 

In both cases we have C2 = 2 [Ho] > . This is especially important in the case 
of quasi-cyclic chains. Consider the solutions of these equations (n = 2) that do 
not depend on y . We obviously have the following solutions 

x = I — = (semi-cyclic case), 

J V C 2 .9 + 4e a / 2 cosh g + C 

x = I ^es-kj W - 9 = /o ' (q uasi - c y cliccase )- 

It is clear that for the corresponding values of the constant C in both cases we 
obtain wide classes of smooth nonsingular periodic solutions. Being independent of 
y , the expressions above determine doubly-periodic operators with arbitrary period 
in y . For the magnetic field we have 

Ho = C 2 - 2e a/2 sinh(/ - a/2) (semi-cyclic case), (31) 

V 2 = 2C 2 -e f ° =e h =H 2 + C 2 , , x 

(quasi-cyclic case). (32) 
H 2 =C 2 -ef° =H, V4 ,y 

To ensure smoothness and nonsingularity of the operator L 2 in the quasi-cyclic 
case, we must have < C 2 . Solutions satisfying this condition certainly exist, as 
can be seen from the formulas. Since the function fo(x) , as well as the potential 
and magnetic field, does not depend on y, it follows that we can choose such a 
(real) gauge so that 

L=(d x +iA 1 ) 2 + (d y +iA 2 ) 2 H + 2V, A 1 =0, A 2 (x)=A 2 , -A' 2 = H{x), 

where L = Lo (semi-cyclic chains) and L — L 2 (quasi-cyclic chains). In both cases, 
after the substitution tp — e lky (p(x, k) , we obtain the following equation in x : 



^kfk = 0, 

A fe = d 2 x - (k + A 2 (x)f + (2V(x) - A' 2 (x)). 



(33) 
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In the Landau case, we have A 2 = Hx , V = H/2 . After that the substitution 
x' = x — kH^ 1 the operator reduces to the harmonic oscillator Ako = d 2 — H 2 x 2 . 
In the case of a nonhomogeneous magnetic field we have H(x) = H 

A 2 = Hx + A°(x), where H= ^ ■ 



TiT 2 ' 



here T\ is the period in x and T 2 is the chosen (arbitrary) period in y , which does 
appear in the formula anyway, because 



[H] = T 2 [ 1 H(x) dx. 
Jo 



The function A^(x) is periodic with period T\ . The substitution x' = x + kH 
reduces this equation to the family of equations 

Afc9?fc = 0, 

A fe = d 2 x , - x' 2 H 2 + U(x' - kW 1 ) - 2kHA° 2 (x), 
U = 2V(x)-H(x)-(A° 2 (x)f. 

The potential V 2 for the operator L 2 in both cases has the form 

V 2 = Vo + C 2 (semi-cyclic case), V 2 = H 2 + C 2 (quasi-cyclic case). 

The quasi-cyclic operators reduce to a one-dimensional family of oscillator-like 
operators: 

A fe = d 2 x , - H 2 x' 2 - 2Hx'A° 2 (x) + U(x), (34) 

where H is a constant and A 2 , U are periodic functions (above). 
Let us consider the doubly periodic nonsingular case. 

Proposition 3. Suppose that the magnetic flux is positive and integer, the point O 
belongs to the spectrum of the operator Lo in £ 2 (R 2 ) and all the operators 11111 
of the semi-cyclic chain are nonsingular, while all the fj, j = 0, 1, . . . , n, are real 
and smooth. Then the point (—C n ) also belongs to the spectrum of the operator L 
in£ 2 (R 2 ). 

The proof of Proposition 3 is based on the fact that the Laplace transformation 
takes magnetic Bloch functions to magnetic Bloch functions. This fact is readily 
verified by substitution and direct calculation. If tp n ^ 0, then the proposition is 
proved. Under its assumptions, if ipn = 0, then it follows that Q n -iipn-i = 0. 
Suppose that j satisfies tpj — and tpj-i ^ 0. We have Qj-itpj-i = 0. Thus the 
operator Lj-i coincides with half of the Pauli operator 

L-=L j - 1 = (E + B J - 1 )(d + A j - 1 ). 

However, we know that the eigenfunctions of the principal state are possible for L_ 
here only for negative magnetic flows [-Hj-i] < 0, whereas by assumption we have 
[Hj-i] = [H ] > 0. This contradiction proves Proposition 3. 
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Proposition 4. Equation (30) of the quasi-cyclic chain of length n = 2 possesses 
smooth solutions, essentially depending on both variables (x, y) , doubly-periodic 

— 1/2 

with respect to the square lattice with any period T > 2irC 2 sufficiently close to 
27rC 2 ~ 1 / 2 . These solutions differ little from f = ln(C2/2) . 

The proof of this proposition (necessary to establish that equation (30) does 
possess a large class of nontrivial doubly periodic solutions) was communicated to 
the authors by S. Kuksin. 

Conjecture. For n ^ 5 , the equation of the quasi-cyclic chain does not have any 
nonconstant periodic solutions close to a constant. The constant is isolated. 

It seems probable that these equations have no solutions already in the linear 
approximation. For n = 2, 3, 4 solutions in the linear approximation can be found. 

Theorem. Suppose we are given a quasi-cyclic chain of length n of operators 
Lq, . . . , L n with smooth nonsingular positive potentials Vq, . . . , V n -\ and nonsin- 
gular magnetic fields H , ■ ■ ■ , H n with positive flow [H ] = • • • = [H n ] , the po- 
tentials and fields being doubly-periodic with period lattice T± , T 2 ■ In this case 
the Schrddinger operator L n possesses two nondegenerate levels 0, +C n each of 
which is isomorphic to the Landau level in a homogeneous magnetic field and zero 
potential. The eigenfunctions of the levels , +C„ can be computed explicitly. 

Proof. Let us construct the eigenfunctions of the levels , +C n . By the quasi- 
cyclicity condition, we have V n — e' n = H n + C n . Hence the operator L = L n — C n 
has the form 

L = L n -C n = (d + A n )(B + B n ). 

Since the flow [H n ] is positive and integer, we can extract the eigenfunctions of the 
ground state from the paper [DN80] 

(L n -C n )i> = 0, V = e v Y[ a ( z ~ • • • > a ( z ~ a m)e az , 

i rr ( 35 ) 

^n = --JJ K ln\a(z-z')\H n (z')d 2 z', z' = (x',y'), z = (x,y). 

Here K is the elementary cell, we do not specify the class a and constants a\, . . . , a m | 
here. The eigenfunctions of the level L n ip = are obtained by Laplace transfor- 
mations from the functions L ipo = . Indeed, we have 

ip n = Qn-i---Qo(ipo), where Q = e^ /2 (<9 + Aj). (36) 

Let us put 

4>o=e ipo f\o-{z-a 3Q )e aoZ , where A<p =-- // In \a(z - z')\ H (z') d 2 z' . 

7=o nJjK 

(37) 

The constant a is related to (oi,. .. ,a m ) by formulas presented in [DN80]. The 
functions ipo, ip\ , . . . , tp n are magnetic Bloch by construction. We already know 
that the level (L n — C n )ip — is isolated, discrete, and isomorphic to the Landau 
level [DN80, N83]. Let us prove this for the level C n . The subspace of functions 
(36) obtained is isomorphic to the Landau level (by construction) and determines 
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a subspacc in the Hilbert space /^(R 2 ) ■ Let us prove that there is nothing more 
at this level. Let us apply the inverse Laplace transformation to the operator L n . 
After n steps we obtain the operator L , while from the functions we get 
ipo . For the level LoV'o = we already know that there is nothing more at this 
level, except the functions found in [DN80]. If the n-fold Laplace transformation 
(ipn) (V'o) is an isomorphism, then the Theorem is proved. If it is not an 
isomorphism, the inverse Laplace transformation has a nontrivial kernel. This 
means, for n > 1 , that for some number j to the operator Lj corresponds a 
function ipj such that Ljipj = and (d + Bj)ipj = 0. Arguing as in the proof of 
Proposition 3, we see that Lj = (d + Aj){d + Bj) . Thus the subchain L 0} . . . ,Lj 
is quasi-cyclic of length j < n, where the constant Cj is zero. Having in mind the 
fact that the magnetic fields are nonsingular and the relation fj = ln Vj , we sec 
that Cj = j [Ho] 7^ ; this contradiction proves the theorem. 

Remark. As the authors were told by Hector de Vega, the equation Ag = 1 — e 9 
was obtained as the "Bogomolmy reduction" of the Ginzburg-Landau equations 
for a special value of their parameter that once divides superconductors of types 
I and II (see an explicit quotation to the work of de Vega and Shaposhnik in 
[NV3]). The applicability of the Ginzburg-Landau equations at this point is not 
clear (possibly it necessitates some corrections, as it was explained to the authors 
by experts from the Landau Institute). However, in this situation the "physical 
magnetic field" differs from ours by a constant (in such a way that the magnetic 
flux in the case of smooth nonsingular solutions is zero). The physical magnetic 
flux in the Ginzburg-Landau case entirely consists of singularities. The study of 
solutions with isolated singularities is meaningful in our case also. In this case 
nontrivial decreasing solutions as \z\ — > oo (z = x + iy) are possible in the case 
C2 = as well, where we have an elliptic Liouville equation, Ag — —e 9 . Its general 
solution is given by the following 

Proposition. For an arbitrary analytic function f(z) the function ip defined as 

- (1 + 1/M1 4 ) 2 { } 

satisfies the following elliptic version of Liouville equation 

Atp + 32e v = 0. (39) 



Proof. It can be checked that the function <p of the form 

eV = (l +r 4)2> 

where r = \z\ 2 satisfies this equation. Now the statement follows from the following 
symmetry of Liouville equation: if ip(z, z) is a solution of (39) then (p(f(z),f(z)) + 
21og|/'(z)| , where f(z) is analytic, is a solution too. 

Taking a rational or elliptic function for f(z) we obtain solutions with singular- 
ities at the poles. The search for solutions with singularities also makes sense in 
the case C2 7^ . 
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Appendix I (S. P. Novikov) 
Difference analogs of Laplace transformations 
and two-dimensional toda lattices 

In the continuous case, the Laplace transformation L — > L , ip — ► was formally 
the same in the elliptic and the hyperbolic case: one could be obtained from the 
other by the formal substitution 

d—>d x , d^dy, 

(B + B)(d + A) + 2V^ (d x + B) (d y + A) + 2V 

(here, of course, the global statements of the problem are entirely different). 

For difference operators, the hyperbolic and selfadjoint elliptic case will already 
differ on the formal level. The separating out of classes of difference operators 
that possess Laplace transformations seems important, since this would indicate 
the most natural classes of operators with maximal "hidden symmetry" (in some 
sense) and therefore with the nicest properties. 

1. Hyperbolic case. Let the operator L be of the form [Kr85] 

Lipn =4> n + a n Tiip n + b n T 2 ip n + c n TiT 2 ip n = 0, (39) 

where n = (n\,n 2 ), Nj e Z, and T\,T 2 are the translations determining the 
lattice, 

TlVVi = Ipn+Ti, T 2 ^ n = ^n+T 2 Tl + T X = (m + 1, n 2 ), Tl + T 2 = (m, U 2 + 1). 

Lemma 5. The operator L has a factorization of the form 

L = /„[(1 + ti„T 1 )(l + v n T 2 ) + w n \. (40) 
The proof of the lemma is elementary. We come to the following formulas 

/„ = (1 + w„) _1 = a n - Tl b n c~^ Ti , u n = u n = a n /f n , v n = b n /f n . (41) 

Note that formulas (41) are local, i.e., to obtain the factorization we need not solve 
any difference equations. 

Definitions. The Laplace transformation is the assignment L — > L, %j) — > ip , where 

14 = 0, Lip — 0, 

~ w„ , ~ (42) 

L= -^[(1 + «„T 2 ) W - 1 (1 + U „T 1 ) + 1], V= (l + v n T 2 )^ n . 
l + w n 

The inverse Laplace transformation is the assignment L L, tp — > ip , where 

Z^=0, i> n = {l + q n T 1 )i> n , L = g n [(l+p n T 2 )(l + q n T 1 ) + s n \. 

Remark 1. Just as in the continuous case, the Laplace transformation and its in- 
verse are indeed inverse to each other. This is verified in an obvious way. 

Remark 2. It is easy to see that to each pair of basis periods (T-j 41 ,!^ 1 ) and 
(T 2 ±1 ,T 1 ±1 ) there corresponds a Laplace transformation quite similar to (40)-(42), 
where we have the pair (T\,T 2 ) and (T 2 ,T\) for the inverse transformation. The 
algebra generated by all Laplace transformations is extremely interesting in this 
case. 

In the gauge equivalence classes L — > gLg^ 1 , the potential w n and the nonphys- 
ical magnetic field H n are invariants, 

w n + 1 = Cn-^a'^b' 1 , e Hn = a n b n + Tx b~^ a~\ T ^. (43) 
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Lemma 6. The Laplace transformation depends only on the invariants of the gauge 
transformation and can be expressed as 

e H " = e H - w n+ T 2 w n+Tl w~ 1 w~l Tl +Ta , 
1 + w n+Tl = e~ Hn (l + w n+T2 ). 

The proof is a direct verification of the formula. 

The infinite chain of Laplace transformations obtained according to formulas 
(44) leads to a completely discrete analog of the two-dimensional Toda lattice for 

(k) (k) 

the magnitudes w„ , where the values of H„ are expressed via (44) 

1 i ( fc ) 

H (k+D _ 1 + W n+ T 2 
1 + W n+T 1 

Here (v)n\Hn^) is obtained by the Laplace transformation from (io^ fc_1 \ Hn 
k G Z, n = (ni,n 2 ) . This system should be compared with [KLWZ]. 

Example 3. Let us consider a cyclic chain of Laplace transformations 

H (k+m) _ rr(k) ,,,(fc+m) _ (fe) 

The case to = 1 is trivial. For m = 2 we obtain a discrete analog of the 
sinh-Gordon equation, and after reduction, we get 

W n +T 1 +T 2 =W~ 1 (C + W n+Tl )(C + W n+T2 )(l + W n+Tl ) _1 (l + I0 n+T2 ) _1 . 

For C — 1 the equation degenerates. 

2. Laplace transformations for selfadjoint real difference operators. The 

ordinary discretization of the Schrodinger operator, say, on a square lattice, where 
at the given point n = (711,712) one takes the (weighted) sum of the values of ip 
at the same point and four its nearest neighbors (m ± T\, 77-2) , (ni, n 2 ± T 2 ) , does 
not possess the factorization necessary for the construction of the Laplace trans- 
formation. The latter can be obtained for the operator L if we use the equilateral 
triangular lattice, where each vertex has 6 nearest neighbors, 

\T 1 ±1 \ = \T 2 ±1 \ = \(T 1 T- 1 ) ±1 l 

The operator L has the form 

L = a n + b n T x + c n T 2 + d^T^ 1 + b n _ Tl T^ + c^T^ 1 + d T2+n T 2 T^\ (45) 

Lemma 7. An operator of the form (45) with real coefficients possesses the follow- 
ing factorization 

L = {x n + y n Ti + z n T 2 )(x n + y n - Tl T^ x + z n - T2 T 2 ~ v ) + w n . 
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The proof is a simple calculation. Note that here, just as in the hyperbolic 
case (40), (41), the factorization is given by local algebraic formulas, which do not 
require the solutions of difference equations. We also have the inverse factorization 

L = (x' n + y' n T^ + z n T^){x' n + y' n _ T J x + z' n _ T2 T 2 ) + w' n . 

To any solution of Lip = a new solution is assigned: 

M = 0, L = (x n + yn-nT^ 1 + y n - T ^ 1 )w- x {x n + y n T x + z n T 2 ) + 1, 

V>n = {X„ + yn-^T^ 1 + Zn-T^ 1 )^. 

We also have an inverse Laplace transformation induced by the inverse factoriza- 
tion. The two transformations are mutually inverse up to a gauge transformation. 

Here the discrete analog of the two-dimensional Toda lattice, as well as those 
of the cyclic, semi-cyclic, and quasicyclic chains arise. They will be studied in a 
forthcoming paper. 

Remark. To each pair of neighboring basis periods of the same length (Ti,T 2 ), 
(T 2 ,T^T 2 ), (T 1 ~ 1 T 2 , Tj -1 ) , (Tf 1 ,^- 1 ), (T 2 -1 ,7iT 2 -1 ) , and (T^Tf 1 ) there 
corresponds one Laplace tranformation. The algebra of these transformations is 
being studied. 

Introducing complex values for the fields, we also obtain the discretization of the 
Schrodinger operator in a magnetic field (what we considered above is a discretiza- 
tion of real operators only, where physical magnetic field is equal to zero.) 

Appendix II (S. P. Novikov, I. A. Taimanov) 
Difference analogs of the harmonic oscillator 

Let us consider the real selfadjoint second order difference operator acting in the 
space C 2 (Z) , n£Z, and given by 

Llp n = V n 1p n + Cn-^n-x + C n 1p n+ i, 

L = c„_iT _1 + c n T + v-l, 

where T: n — > n + 1 is the translation operator along the lattice. Here c n ,v n e K, 
T+ =T-\ and L+ =L. 

As in the continuous case, one defines the factorization of the operator L : 

(I) L + a={a n + b n _ l T- 1 ){a n + b n T) or 

(II) L + a = ( Pn + q n T)(p n + ^-iT" 1 ). 

Even without requiring the coefficients to be real, such a factorization will exist 
for any a . In order to find the values of a n , b n , p n , and q n from v n and c„ , 
we arrive at the difference analog of the Riccati equation, which appears when we 
factorize the continuous Schrodinger operator of the form 



—d% + u(x) + a = 



-(d x + v)(d x -v), v 2 + v x = a + u(x). 



EXACTLY SOLVABLE TWO-DIMENSIONAL SCHRODINGER OPERATORS 19 
In the discrete case we have 

(I) v n + a = a 2 n + b\_ x , c n = a n b n , 

(II) v n + a = p 2 n + q*_ u c n = q n p n . 

If all the coefficients a n and b n are real, then 

(I) L + a = QQ+ 1 Q+=a n + b n T. 
Similarly, if all the coefficients p n and q n are real, 

(II) L + a = RR+, R+=p n + q n T-\ 

The Darboux-Backlund transformations B a , by definition, are given by 

(I) L -> Q+Q = 1 = B^L, 

(II) L -» R+R = L = B { l T) L. 

It is obvious that these transformations can be taken to be inverse to each other 
after an appropriate factorization: 

Bf\B^L) = L + a, Bi I \B^)L)=L + a. 

(Recall that in contrast with the two-dimensional case, the factorization here re- 
quires solving a Riccati type equation and is therefore nonunique.) The cyclic 
chains, just as in the continuous case, are determined from the condition: 

Ln — B aN ■ ■ ■ B ao L, Lj = B aj Lj_i, L = L^. 

It is not difficult, by following [VS1], to prove the following 

Proposition 1. Cyclic chains satisfying the condition X^Li a j = consist of 
finite-zone difference operators Lj corresponding to one and the same Riemann 
surface of genus g < [N/2] (see [DMN, DKN2]). 

More interesting are the chains for which ctj = h> 0. 

The first difference analog of the harmonic oscillator. Already for N = 1 
nontrivial phenomena arise. Consider the cyclicity conditions 

L = QQ+-h 1 L = Q+Q = B h L, Q+Q = QQ+ + h, Q+=a n + b n T. (1) 

From (1) we obtain an equation of Riccati type, which implies 

a n = a = const, b n — b n _ l = h (2) 

or b n = \/nh + b1. Suppose that further 

o = l. (3) 

We come to the following conclusion: the operator L = QQ + is not defined in 
Li(/£) as a real operator. 
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Proposition 2. 1 The operator L = QQ + described by formulas (l)-(3) deter- 
mines a real self adjoint operator L in the space Hi C £ 2 (Z) if and only if I £ Z 
(the quantization condition). 

(Here ip n G Hi if and only if ip n = for n ^ —by/h = I.) The space Hi is 
isomorphic to £ 2 (Z+) with zero boundary condition for k = 0. 

Proof. If the number I is an integer, then 



Qt = (l±y/(n-l)hT). 

Setting n — I + k , we come to Z + . It is readily verified by substitution that the 
operators L± = Q±Q± take the space Hi to itself if and only if I is an integer. 
We come to the following operators in £ 2 (Z + ) : 

L± = Q±Q±, Ql = (l±VhkT). 

The spectrum of the operator L± = Q±Q± in the space £ 2 (Z+) is the same as 
that of the ordinary harmonic oscillator. The ground state will be 



L±*l>± = 0, Q^o = 0, V? fc 



(±l) k - 1 /y/h k - 1 (k-l)\, k^l, 

0, k < 0. 



The higher eigcnfunctions have the form 

Lemma 1. The eigenf unctions tp^ of the operator L± have the form 

.ffere i/ie Pj(fc) are polynomials such that 

P j (k) = (l-h(k-l)T- 1 )P j _ 1 (k), P = l, and B(k) = j ' 



1, fc>l, 
fc < 0. 



The expression is a Poisson distribution, while the polynomials Pj (k) are 
orthogonal with Poisson weight tp^(k) on Z+ . Undoubtedly they are known, 
although their relationship with harmonic oscillators and bosonic commutation re- 
lations, most probably, was never discussed. 2 



1 Similar operators were constructed in [AS90]. However, the analysis of the corresponding 
formulas does not reveal if their authors have the same thing in mind as we do. Apparently our 
results are new, at least methodically. 

2 Thcse polynomials arc known as the Charlet polynomials. 
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The second difference analog of the harmonic oscillator. Let N = 1 . Put 

L = QQ + and L = Q + Q. We shall introduce a family of operators L depending 
on two constants c, a e R : 

L(c,a) = Q(c,a)Q + (c,a), L(c,a) = Q + (c,a)Q(c,a) 

so as to have 

a 2 Q+(c,a)Q(c,a) = Q(ca 2 ,a)Q+(ca 2 ,a) + D, D = a 2 -l. (4) 
Let us put 

Q+ = l + ca n T, Q = 1 + ca n - 1 T-\ a ^ 0, c ^ 0. 

Consider the transformation r : n — > 1 — n . We have the formula 

TQ(c,a) =Q+(c,a- 1 )T. (5) 

The relations (4)-(5) are extremely interesting. Apparently, they have not appeared 
previously. 

Theorem. The spectrum of the operator L(c,a) for A in the semi-interval [0,1) 
in the Hilbert space £ 2 (Z) has the form 

1) a > 1, A„ = 1 - a" 2 ™, n > 0, 

2) a < 1, A„ = 1 - a 2 ™, n > 1. 

TTie eigenf unctions are the following: 

1) a > 1, Vofe(c, a) = {-l) k c- k a-( k -V k / 2 , k e Z, 

V»n(c, a) = Q(c, a)Q(ca 2 , o) • • • Q(ca 2 "™ 2 , a) V>o(ca 2 ™, a), 
V'nfc (c, a) = P n (fc, c, a) ip k (c, a) , 

P n (fc, c, a) = (a- 2fc - ca 2 "- 2 )™, n > 1, 

2) a<l, Vi(c,a) = Tip (-^, -), 

^n(c,o) = Q+(^,a) •• -Q + (^2,a) ^1(^=2, a), »» > 2. 

Remarks. 1. For A ^ 1 the spectrum of the operators L(c,a) is not known to the 
authors. We conjecture that it is continuous. 

2. If a > 1 and c < , then the operator L may be considered in £2^) so that 
its restriction to the family of lattices i + Zcl yields our family L(c, a) : 

J L_(m, 7 ) = (l-7 a:_m_1 P _1 )(l-7 x_m P), 7>0, meZ. 
The ground state ^0(2) {L&o(%) = 0) is a function satisfying 

(1 - j x ~ m T)<P (x) = or <P Q (x) =j x - m <P (x + l). 



22 



S. P. NOVIKOV AND A. P. VESELOV 



It has the form 

$ (x) = 9(x)'y x{2m+1 - x)/2 , 

where g{x) = g{x + 1) is any function of period 1 . 

Let us normalize the ground state by defining g{x) from a continuum of normaliz- 
ing conditions so that the expression <£q (x) is a normalized probability distribution 
on any lattice of the form 6 + Z C M : 



i ,'ax{m + 1) ax 2 x 
g (x) = exp ( — ] 6 



a(m + 1) 
ax 



where a = 2 In 7 > and Q[u | a] is a theta function: 

S[u\a] = J] cx p (- ^j- + 



nu . 



It follows from the analytical properties of theta functions ([BE]) that g{x) is a 
smooth function with period 1 . The normed ground state is given by the formula 

x(x-l)/2 

*o(s) 



Q[(m+ 1 - 2x) In 7 | 2 In 7] ' 
Thus the cyclic chains are of two types: 

1) L = B aN ■ ■ ■ B ai L = L' , where L acts in £,2{1> + ) for an appropriate "quan- 
tization" of the parameters; 

2) the operator L' coincides with L after multiplication by a constant and a 
translation along x in the natural realization in £2 (M) , where the restriction to 
the lattices {5 + 1) c K generates a family of discrete operators participating in 
the definition of the cyclic chain similarly to the case N = 1 considered above. 

Added in proof. In the paper [SVZ], the problem was essentially posed already. 
However, that paper contains unmotivated restrictions. For example, the constant 
5 = In I c|/ In a is assumed rational in [SVZ]. Further, the assertion in [SVZ] accord- 
ing to which the relation (4) above or the corresponding relation (14) in [SVZ] and 
its consequences (15)— (17) "clearly define a spectrum generating algebra" is incor- 
rect. This assertion is certainly false, for the case in which tp satisfies the equation 
Q + ipo — and growth exponentially. This is indeed the case when q > 1 , a situ- 
ation omitted in the papers [AS91, SVZ]. In this situation the authors have found 
a different ground state with eigenvalue 1 — q {q^ 1 = a 2 ) by using the additional 
symmetry r (see Theorem 2). 

Moreover, the assertion that the relations mentioned above "clearly define a 
spectrum generating algebra" is incorrect for another reason: it gives no information 
on the spectrum of the Schrodinger operator L = QQ + with A ^ 1 in the Hilbcrt 
space £2(2). According to our conjecture (see Remark 1 above) this spectrum is 
continuous and occupies the entire strip A ^ 1 . 
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